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fUPULA HON PROJECTION * 

1 . IMRODUC1 ION 

Vshat will be the distribution of ages in the U.S. 
population lu >ears< from now, SO years from now, or' 100 

from now' 7 Many things depend on the answer 'to that 
quest 1 on- - the future of the social security system, the 
demand for health care ser\i«ces, the enrollment in colleges 
and un iv ers i 1 1 e s to name a few If we know the present 
age distribution of a population i sa> , how many individ- 
uals are •■-+ years of age, 5-9" years of age, etc.j can we 
predict the age distribution in 50 years 0 It seems clear 
th^t an> such -pred ict ion would be subject to lots of 
^rror. \ftvr all, we cannot foresee future change^ in ^ 
birth and death ra»tes, nor ^an we estimate the effects of 
migration. \ 

\ somewhat simpler question which can be easily for- 
mulated and answered mathematically is this: If current 
birth rates and survival rates continue unchanged, what 
will, tne population distribution be at a future 'date 7 (assum 
i ng the effects of migration are ignored) 0 .Although this 
question is simpler, it is not withcut practical interest. 
J or example, it is of value to explore mathematically the 
conseuuences of various birth and death rato> without hav- 
ing to wait until tne popula t i on "ac tua 1 1 y experiences these 
'con?cqucnet^ . " * ' 

Several mathematical model* have been developed for 
thi- simpler question. Here we shall examine one of these - 
jl mitrix model devised hv F ; \\. Leslie and others^ In this 
m/del, matrix multiplication i^ used to update tne popula- 
tion from one time^period to another. As we shall see, 
tht stuJ, c 1 such a pro^c^s leads us to examine powers of 
mat rice- hir object will be to learn something oi the 
Scha\ i r r of p (S ^ers oi matrices and to sqc how knowledge 



of eigenval ues *and eigenvectors helps in understanding 
t his behav id*r . 

2. THE MODEL 

2.1 population Proj ec 1 1 on - - An Example 

Let's begin a simple example. Con-s^rtier a'fic- 

titious animal population consisting of 100U ydung animals 
l-l vedr old, 8« n individuals 1-2 years old, and bOO indi- 
viduals 2*-3 years old. \se will assume that none of these 
aniTi-als lives longer, than 5 years, «He can record this 

7l0un^ 
I 800 1. ; 

I ooo. 



\ or rrLvitv *>e shall cal'l the ate groups Class I ■ 0-1 year J , 

At 



Class II . 1 - 2 years; , and Class fl I i2-3 vcarsf. % 

. f I 

t j" 

1 ve wan^' to know the age^ d l s t r i but i on vector one vcar 
iron 1 now we will need to .know two things 

* * 

1/ the proportion of thos'e animals currently ajiv^ 
n> ( " that ■ wi 11 -survive until next year, and 

^5 2. 'how man) offspring will he born and will survive 

long enough to* r>e counted next vcar. 



In our example, I^e-t^u, «. uprose that 1/2 of^the individuals 
in Class I d^c , 5 U f ' individuals) survive to be in Class » 
II the next year, and let us also suppose that 1/2 of the 
individuals -in Clas*s II f i , e . , -4 00 individuals} survive to 
be; in CI as*. III the next v;ear. f Under our assumption, 'indi- 
viduals currently in. Class Ill-will be dead by ne*t year. J 
I his process is indicated by the solid lines m Table 1, 
In addition let us suppose that individuals in Class I pio- 
dxi l 1 0 no offspring, that ea'ch individual in Class II produces 
offspring on the average, and that each individual m Class 
•III produces 2 offspring' on the average; (Here wc are 
including onlv those offspring who survive long enough to 



TAB LB 1 



Class 


NOW 

Number of tndivixiuals 


NfiCT YEAR 
Number of Individuals 


Class 


I 


1000 - — 


r ^l(800) * 2(600} = 2000 


I 


II 




4(1000) = 500 


II 


III 


600-""" ~ 


"-^4(800} = 400 


III 



be counted the next year.} In Class I the next >ear vse 
uould, expect to ^have 800 young (produced by last year's 
Class II individuals) and 2(600) = 1200 young (produced 
by last year's Class III individuals) ,< for a total of - 
2000 nev. Class I individuals. This process is indicated^ 
.by the dotted 'lines in Table 1. r 

T*he .computa t ions we have done can be summarized in 

this matrix computation: 

„ % 

'2000] % 

! soo' 



(0 1 2^ flOOOl 

J 1 ! 

! H 0 0 ' i 800 , 



(0 \ Oj { 600 j 



400 J 



This computation is of the form Ax 



Q 



vshere A is 



the mat r i-x^j:tfnt a in ing the birth and survival paraflfete r s , * 
Xq is the initial age distribution vector, and Xj is the 
age distribution vector after*one year. We u i 1 1 ca41 the ' 
matrix A a Leslie \atrix. * ^ 

2.2 Extending the Projection 

By ,per f ormmg a matrix multiplication v*c have found 
the age distribution after one year. What if v*e uanted to 
knou the age distribution after 2 years 0 If v*e believe 
that birth and survival ratcs'will remain unchanged, then* 
v-e can again multiply by the Leslie matrix A, this time 
using x^ as our starting distribution. Thus v*e compute 

'ft) 1 21 f 2000 f ( fl 300l 
\\ 0 0; I 500 L = I l(M)o! = x; . 
(0 h 0j i 400J I 2 SO | 
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Notice that x 7 ? Ax ^ = A(Ax Q j = A x Q , 



We needn't stop with 



vears of course: 
JL . 



He could compute, 



A^x f 



etc . 



^x^ = Ax , = \ ( A x q j t /\ aq, 

More generally, if the birth ancr* survival parameters remain 
constant lor if wish to know what w^uld happen if they 
u.t remain constant], v,e could find thj^* age distribution"* 
sector after k years {let's call it by computing x^, = 

K V 



hxerc ij>e 1 . Continuing the example of this section^^ojnpute x^, x^ , 
x. and x^. (Note You will probably find it easier to obtain x^ by 
computing Ax,, rather than first computing \* and then computing A^Xq.J 
Do >ou notice aay qualitative trends (e.g., does one class have con- 
sistently rwrc or fewer indijviduata than others}? Try the process 
again using the same Leslie matrixVut a different starting vector'x . 
Does the choice of starting vector seem to affect the trend'' 



v 1 -S\ . ) l SL F_^_ llL Jj 1 c lJ ^- 3 ic Matrix 

The exercise raises 'tuo Questions: 

l- e Jfthe multiplicative process is repeated again 

and^gain does the distribution of ages change ran- 
domly or is there sonic recogni zable* pat t ern *in the 
* successive age distribution vectors? 

Z. Dees the «rttimate benavior of the age distribution 
vector depend on the initial distribution 9 

If there arc any patterns they should show up as l v%e ' 
compute higher and higher powers of the Leslie matrix. Let's 
cqntmue our exampl'e'by computing some powers of A. \ou can 
do this by hand fa bit ted'ious, of coursek or by computer (a 
program for d6mg this is shovJjpin Appendi^c A). Here are 
some results (numbers are rounded to 6 decimalsj: 



A = 



( 0 


1 21 






'.5 


1 




. 5 


0 0 


A* ■ 




0 


. 5 


1 


0 


.5 0, 






.25 


0 
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,4 



In 



'.25 .1 I j 
.25 .25 .5 
..125.25 l)j 

' .402344 .:y'6S75 
' . 199219* .402544 

. .'199009 .199219 



.4 375 
{ .1875 
^.0937 



*75 75 1 
.4375 .375 



. 7 9 08 " S | 
.39843*' 
.2031251 

A*' J 



o .1875 .25 j 
f. 40(^09°' .799988 .799988 



j . 199997 
( .099999 



.400009 .399995 
.199997 .*00012 



4 64 




I . 4UU00C . SnOOOH . 800000] 
| r 200000 . 4H0U0U .400nnj 
i . l nnooo . 2uoouo . 200000 )% 

We sec that powers of A are much alike for "large" powers 
'at least this is true in our example--we will come to a more 
general case later j. V»hat does th'is observation impfy about 




TD age distribution vector'' To see, let's suppose the 
initial age distribution vector is 



^0 



l b ! 



Then if V 

J v 



is lai-ge, x, 
/ 

{ .4 .8 . 8Ha1, 
I i 1 

. .2 .4 .4Mb 



( =A x n ) is approximately 



where x 



I 

t .1 
= a - 



2b 



2jicj 
+ 2c. 



v 0 

f.4a + 

..a 
[!la * 



.8b 
.4b 
. 2b 



.8c) 
.4cj 
.2cj 



. 4 I 



The results are rather surprising* Although the total 
size of the population depends on the initial values a, b 
and c, the relative proportions of individuals^ in the three 
age classes approach fixed, rAios 4:2:1 as k+*>, and these 
proportions do not depend on the initial age distribution 
of the population. (This is sometimes described by saying 
that the population "forgets" its initial age structure.) 

Although the ratios 4*2:1 are never quite reached in 
a finite length of time, it is interest ing *o notice that 
if th£ ~}>opul at ion did achieve these proportions at some time 
m, then the age distribution vector would not change during 
subsequent time periods. Tor example, if 



9 

ERIC 



'J 



40 
20 
[10J 



then 





f0 1 2' 




'40' 




•40 1 


A % = 


!j 0 0 




to 




20 


,0 4 0. 




.10; 




. 10 . 



1 



• The, fact that A.x m = x can be expressed mathematically 

mm r - 

by sa>'*r7lg that x m is an eigenvector corresponding to the 

eigenvalue 1. (Recall that a* nonzero column vector v is 

an e igenveozor of a square matrix A Tf there is k scalar 

>, called an eigenvalue > such that Av = >v.J \ 

In case it escaped your notice, go back and observe 

that the' ratias 4:2:1 which occur m the eigenvector also 

•> occur approximately in the columns of A when k is a large 
64 

number 1 see A ^ for instance) . 

2.4 Another Example * . 

Based on the example of the preceding section v^e have 

several hunc hes * abo ut possible theorems. Before exploring 

these hunches it vsould be wise to look at one more example. 

si 

he need to knovv whether the behavior shown^in the previous 
example was typical of population growth using Leslie '* 
matrices. Consider the Leslie matrix 



B 



1 4 

h 0 



fCan you interpret the entries of this matrix using "popu- 
1 

latioft" 1 a.nguage ? ) ^ Some powers of B are shown: 



B 



B • = 



5 4 
.5 2 J 

171 340 
42.5 86 



B* 



,16 



11 20 
2.5 6j 

43691 87380 
10922.5 21846 



Here v*e have a real, population explosion! At first glance, 
the clear-cut patterns observed in our earlier example seem 



to be missing. But if you look closely at B and B 



16 



you 



\ 



will see that the ratio of each first-row number to the 
corresponding second-row number is* approximately 4' to 1. 
Taking our cue from the example of Section 2.3, we might 
wonder" whether a vector having a 4 to 1 ratio is an eigen- 
vector of B. By 'computing 



fl 4 

{h o 



41 

I 

1J 



the answer is yes 



° = 2, 
2j - [1 

However 



here the eigen- 



v*e see 

value is k, and that- gives a clue to the' reason for the 
population explosion* even if the ratio in the age dis- 
tribution vector stabilized at 4 to 1, the population 
would continue to grow and would, in fact, double every 
time period. 

The. analysis of the powers of B might have been 
clearer had we accounted for the doubling tendency of the 



population by dividing e*ntries of B 
and so on. Fibre are-the results: 

r 



by 2 



2 



of B 4 by 2 4 : 



1 R 2 



1 



.125 



r 



1 D 16 



.6875 1 . 25 J 
.15625 . 375 



,16* 



. 66667 2 1 . 33331 3 
.166664 * \ 333344 



1 D 8 _ f. 667969 1.328125 

2 [. 166016 .3*35938 

(Computations are rounded to six decimals.) The patter^ 
now is much , clearer , and oirr computation v*ould lead us to 
guess that* „ * 



lim i^B 



-?/3 4/3 
1/6 1/3 



f 



Of course, computation, however useful, is no substitute 

for understanding. How can we analyze the behavior of A ? 

Can we predict the ultimate form of A without actually*. 

computing the powers? * ~* 

* { % <_ ,. . 

By lim C we meaiuthe matrix, if it % exists, whose (i,j)-th entry is 
k-*<» ' » 

. (k) * f , . y 

, c iJ , where c\J denotes the (i,j)-th entry of C . 



C 



Our observations up to this point suggest that a knowl- 
edge of eigenvalues and, eigenvectors would be valuable in 
oigr analysis. We already know one eigenvalue and a cor- 
responding eigenvector for the matrix B. Let's find the* 
other eigenvalue/ Recall that the eigenvalue's of B will 
be roo.ts of the characteristic polynomial' 



f 

det(H - B) =■ det| \ 
I 



1 0 
0 1 



fl 4 

0 



det 



= *>** - A - 2 - (A - 2) ( A + 1) . 

Ffrere are tv,b eigenvalues, \^ = 2 and A-, = -1. An eigen- 
Vj corresponding to eigenvalue Aj can be found by 



vector 

solving the equation (A 
matrix equation is 

',1 - I 1 'X | fol 

i . ! = ; i 

One so 1 ut l on for v 



j- B)v j 



0. For A 



this 



- >>- 4y = 0 
^x + 2y = Q 4 



where v 



1 



1 



i: 



CAn> multiple 
. responding to 



eigenvector corresponding' to > 
> in liar manner , v 



f this vector is also an eigenvector cor - 
j--the Gig ens pace i^ one -d lmensional . ) An 
-1 can be found \n a 



i-ij 



and 



0 



Ihe general theory of eigenvalue^ tells us that if we let 

v, and be the columns o-f a matrix P, then P u i 1 1 be 

j • - 1 

in\ertible and POP = B ^ where - » / 

p = ll-lj 

Of what use is ?his ;n computing powers of B ? We 
see that 

'1 _ ™2 n -l 



1 

0 > 



0 -1 



PDP 



PD~P 



B 5 = P^P* 1 



B k « PD k P* 1 



This observation is useful because *\power^ of a diagonal 
matrix are very easy to compute, ^n oui***example 



{-!)' 



It follows that 



4 2 
1 -1 



1 -1 



2* 0 * 

0 (-1)£ 

2 k D 

o (,i) k 



4 2 
1 -1 

1 . 1 
2 



-1 



It is helpful to rewrite this in^ the foi 



1 R k . 
2 K 



4 2. 



1 -1 



1 0 



i a 

1 2 
6 3 



Since lim = 0, v.e see that 

A ... 

4 2 



1 lm - 
k-*°° 2 



v* 1 



i -i 



1 0 

0 0 



V 2 

which checks with our initial computations.- 



1 1 
S 3 



2 4 
1 1 



Now that we have determined t i:he limiting matrix for 
k k 

(1* )B , whatv can we say about the age distribution vec- 

tor^ ff^ . 



is the m'itial age distribution vector, then )x, , which 
k * k 

equal ^(^ )(B x^), approaches v * 
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J 



fa 4] 
3 3 

i i 

U 3j 



= (a + 2b) 



2 
I 

1 
5 



4/6 
1/6 



when k 



as k+». Thus x k is proximately 2 (a + 2b) # 

is large."* Notice agairi that ?he. relative proportions in 
the, age classes are eventually almost const-ant (801 in 
Class I, 201 in Class II--ratio of 4 to 1). Also notice 
that this 4 to 1 ratio is" found in' the eigenvector cor- * 
responding to Aj = 2. The di f ference^ between this example^ 
and the example of Section 2.3 is that now the total popu- 
lation grous by % factor of approximately 2 each, time k 
increases by 1. 



Exercise 



Let 



!l 3 



A = 



a. Interpret the entries of A in terms of births and survivals. 

b'. Compute the eigenvalues and eigenvectors of A. 
k k 

c. Determine lim A /\ , where \^ is the larger of the eigenvalues^ 



of A. 

d. Approximately what percent of the population is in Class I an<^ 
^ whatipercent in Class II eventually? 

e. How fast does the total population grow each time period when k 
i$ large 9 > 



T 



THEORETICAL BACKGROUND 



3.y Some Observations 

In all the examples and exercises of the preceding 
sections there were certain common characteristics: ' 

V 1. The age distribution vector eventually bphaved 
like a multiple of some fixed vector. 



1/) ' 
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2. This-fi;ced vector was an eigenvector corresponding 

•to the largest eigenvalue. 
'3. 'This largest eigenvalue was real and positive. 
4. The population eventually tended to grow at a rate 
equal to this largest eigenvalue. 

The natural question now 1s: "How t>pical were these exam- 
ples?" Could vve expect this behavior to hoid in general 
(probably too much m tg^ expec t ) or under what conditions u:ald 
it 'hold? . - 

3.2 The Perron-Frobenius Theorem 

Since everything seems to de*pend on having a positive 
--eigenvalue which is'^larger. than the other eigenvalues, our 
first question might logically be: "Is there always such 
an 'c igen\ alue°" Also, can we fi-nd a nonfiega t i ve w *eigenvec- 
tor corresponding to this eigenvalue? (After all, to 
realistic, age distribution vectors must have nonnegat 
components. j It is' surprising that with only some mil 
assumptions the answer to these questions is "yes." A 
famous theorem, known as the Perron-Frobenius Theorem goes 
a long way Toward answering these question?. Jhe proof of 
this theorem is4>e;ond w^at we intend to do here, J but we 
shall at least see how the theorem applies to Leslie matrices. 

Recall that a Leslie matrix has the form 

f a l a 2 a 3 a n-l a n] * 

j bj 0 o ... 0 0 \ 



A = (0 b 2 0 . . . t) ^ I (I) 
1 0 0 b, ; . . 0 0 ' „ 

in o o .*.-. b n .J < j 



"\ KTnKtj + iiJi, sector is one for which all components are nonnegat lve. If 
all components are positive, the vector is called a pocit'jt. vector. 

^If you want td see more abq.ut this theorem you can consult the book by 
Varga fsee the references). 



1 1 



ivhere a j fl ~ •*•» n * gives vt he average -number of off 4 
spring born to an individual in the ith class and 
b^ii =1, n-lj gives the probability that an indi- 

vidual in the ith cl^ss will'sufVive to be counted in the 
f i + I ) > t ^las< during the next tijpe period. Naturally we 
assune that a ^ 0 for l s 1, .... n.' We m 1.1 assume that 
the'survival probabilities are po|itive, i.e., > 0 
v i = 1, . . . H n- 1 ) . Io satisfy tne hypotheses of the 
Perron-f robemus Theorem, we must ^also assume that is 
strict 1> positive, This^ assures tnat the Leslie matrix 
i > lrrvJuc lbje- -a concept v>e shall not explore here. Under 
the^e assumptions, £he Perrdn-Fr^obenius Theorem guarantees 
tne following : w 

1 . hh^ma^nx -V has a posrti^e eigenvalue^ call it 

j, such that ~ '' x \ fp***all other eigenvalues 
. • i» can be real cV complexj. 
^ corresponding to -there exists an eigenvector, 
call it having aj 1 positive components. * 
>. Jhe eigenspace co rjr^ej>n>Dnd l ng to- >^ is one- 

Uinensional , i*.e. , any eigenvector corresponding 
to '| is a multiple of'v^. 



fxcrvise > fhe structure of a Leslie matrix is rather simple and 
\ 

•makes it possible to make some conclusfons about ei-genvalues and 
eigenvectors without appealing to the jierron-Frohenius Theorem 

a Write out the characteristic equation for an arbitrary 3 ^ 3 

Leslie matrix of the farm (Lj , Using Descartes' rule, of signs, 4 what 

can you sa> about the number of positive eigenvalues 9 What about an 

aibitrary n x n Leslie matrix 9 

fx 1 * 
b If yi ls^^o be an eigenvector of Leslie matrix A corre- 



sponding to eigenvalue ^, it must satisfy 

hi ° ° iH = 

* • 10 b 7 0 J(zJ \ 



4 IS you are not familiar with Descartes' rule of signs, you may loflk 
it up in Uspensky,* J V. , Theor y frf bqua 1 1 on s , McGraw-Hill, 1948 

^ 12 



^ " r 



/ 



Using the second t»o equations s oHe for y and z .in terms of x (x ma> 
be assigned an> value, say x = 1], and in this way obtain a fornix 
for a positive eigenvector correspond ing to hxplain why the 

values of x, > and : whicb >>ou found automatically satisfv the 
remaining equation. a } x * a^y * a,: = * jX 9 Can you write a formula 
► for a positive eigenvector in the n ' n case 7 



— ~ ~ " 7 J * 

3. > \rf hx'jmple with Oscillations • 

• , -\s you can see, tfre Perron-pjrobcnius Theorem gives us 

a:^<?<j; what Ike want However, in Analyzing the examples 

of Section^ -.3 and 2,4 it was important that for each 

i » /* 3 ) k * f> as Thlb required, that ' he 

c i i 1 i 

strictlv le>s than 1, i.e., that >, be ztvijil. greater 
f 

zhah the modulus of any other eigenvalue. fWhen *j 1 / 
^r all 6ther eigenvalues / we will say that is. 

io^lKjint . y Unfortunately, the Pe r ron %F ro_ben*ius 
Theorem guarantees only that ^ ^ J >* } ! , and^in fact & 
little experimentation shows that £he assumptions we* made 
fo*r that theorem ax;e not sufficient to justifv strict 
inequalit). I r . tins exercise 



I xemse 4 . Consider the Leslie matrix 

* To si 

A = » J ■ 

rioi 

Starting with the initial age distribution vector x Q = j i0 j, compute 
^x* , X of x„, . , x & . What do you^^erve 7 Find the eigenvalues and 
eigenvectors of \, Then try .to out the limiting argument used 

in Section 2 . J . Mat happens 0 



In tire exercise we found 'that ! ' j 1" ^ : A flc ' somet hi nj; 
hT-xh d i^tl not happen in our previous examples. A quick 
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1 ^ 

X 4 



review of the methods used in our other examples shows 
that the limiting argument that we used breaks down for 
the example in bxercise 4. As we see, the existence of a 
second eigenvalue whose modulus equals >j leads to oscil- 
lating behawor in the age distribution vectors 1 . 

There are rather simple condition^ which will guaran- 
tee t\r*£ is strictl> dominant. One such condition is 
used in* t i is thtcie-m A 

If 'i Leslie ratrn of t^e for" 4 iLj aatisfie> the 

conditions of Section, 3. 2 and if in addition 

there are two consecutive indices i* and l * 1 such 
— r *— 

that a and a , arc K otn p&Mti \ e, then is 
1 i*l — - -- ■ 1 — ' i 

> t n c 1 1 ; . dominant . 5 

( -r f « r ^ of populations, t h i s t h ec r err. requires c H a t there 
r>e i a o v. i^-.sUt .vc a ^ e cla?st> raving p c 1 1 : ve fertility--* 
j r^'iirft.'.t t Vi t A^ll as ua 1 1 * he met in practice. 

h-'.? iii^ a^ain of tne c.«a~ples of Section 2.5 and 2.4, 
enc final question ^i fc K t arise. In thc>e example? *e *ere 
trie t > d i agon a 1 1 ft ' t he I c ^ 1 1 e Tutrix, h } a t if « ecu 1 dn ' t 
ua^a; I c rtunatel . , tne conclusion- about li-ut.r.i; 



a; lor jf f.e arc distribution- \ectoi do • net ncid or 



a : DC in*; i K 1 e t j dugcalizc 1 1 1 e Matrix. I i . o i ire 
far i liar Ait." the Jordan wJ n or.ual forr% ; on van ->ee rr. « 
to extend our ir.al.sis to tnc more general east 'At will 
not fur sue t* e nond l age na li zabl e ^a.^c^ !Crc - * 
3.4 ^UTirv z 



/It'- Aortrwhil*- to stop ana pull t o^ e t n e r trie t hoc r » 
hi.n w «. h i* \ t» o ? i r . first let'- r e c a 1 1 * o u r as- - A rri^ r ion- . 

\ V nT.j>t j_p _n> * h v 'po pu 1 a 1 1 on is go v < r ned b v a I e s 1 i e 
"tu t r l a c i the f o i n i L j i n which 



'if you ire interested in the proof of this theorem, see the hook by 
Pol i lrJ nentioncd in the references \ fascinating (and more central 
tneorum require^ only that the greKc^t tomnon divisor oi the set ot 
lndlvCb i for whi^h a^ - 0 be Qjftc 



I i 



1 . a, , a 7 ^ 



a n _j are nonnegative 



a and b, 
n 1 



b , are positive, 
n - 1 r 



and 



at least two consecutively indexed a-'s are 
po i it ive . 

Hopefully you are convinced bv now jjalthough we have 
given no formal proof) that the -method of Section 2.4 can 
be carried out for any matrix sat isf yin&^pur assumptions 
and that the following conclusions will nold: 

Cone lus ions : As k--°, the proportions of mdi^vi- 
vi duals in the various age classes of approach 
fixed values, and these values are determined by 
a posit ive eigenvector corresponding t (jit he domi- 
nant eigenvalue. The dominant eigenvalue gives 
the eventual growth rate of the population. 

5.5 A simplification Using Left and Right Eigenvectors 



v in the examples we v saw that A equals 
«'l 0 'J ... 0 ^ 



i f ' 2 
n v 



e 



f,. 



\ 



i K ; 

h 4 



which for k large, is approximately 

p l l °] p" 1 

i0 Oj 



where O's represent blocks zeros of the appropriate 
size 1 . In other words, the effect of the eigenvalues other 



than the dominant" one , 



' 1 1 



is less and less as k grows 



larger. V»e might wonder whether we could approximate A 
without paving to compute the smaller eigenvalues for their 
eigenvectors; at all. It tu*Vns out that we car, do this. 
We will state" the facts that we need here. If >ou want to 
understand the justification for these facts, see Appendix 
B . 15 
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In order to <>tfate the desired conclusions we must 
fir.st mention the eigenvalues and eigenvectors of A 7 , the 
transpose of the Leslie matrix A. Perhaps you^already 
know- from your study of linear algebra that A and A 7 have 
the same eigen\ alues . «*I f not, do this exercise. 

Exercise .5 . ket A be an n x*n matrix. Prove that A and A . have the 
same characterise ic polynomial and hence the same eigenvalues. Give 
an example to show that A and A do not in general have the same eigen- 



Let's turn again to our analysis of powers of a Leslie 

matrix A. Suppose ^ is the dominant eigenvalue of A. Of 

course, b> the exercise above* is also the dominant* 
j i » 

eigenvalue of A . Let v De a positive eigenvector of A 
corresponding to and let u, be a>positive eigenvector 

r 6 

cf A corresponding to ^. Th* theorem we have is this 

Under the same assumptions stated m Section 3,4, 

■ k k T 1 

A tends to v ^u-^/u^v^ as 4 ]<-»*> . 

It is important to examine the dimensions of the products 

in v 'l u r /u l v l" *f ^ is p x n/then both Uj asd are n x 1 

vectors. So u{v is 1 x 1 and so can be treated as*a 

<calar. However*, v^u^ is the product of an n x 1 matrix 

b> a 1 fx n matrix and so has dimensions n x n (the same as J 

\, as it shciild; . Incidentally, you perhaps noticed in 

k k f 

previous examples that A //, tended to a matrix of rank one. 
Can you convince yourseli that VjUj always has rank one? 

Now^what about x^? We see that for k large, 

6~~ " ~ " " 

uj is sonetimes called a „eft eigenvector of A. The reason, by 

definition^ u] satisfies A T Uj = AjU^ but this equation is equivalent 
to u^X = 'luf, hence the reference to left eigenvector. As you might 
guess, is sometimes referred to as a right eigenvector of A. 
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k Xu >i f„{ Vl j °. 1 ; 

i 

Since UjX,, is also lxl, you can see that is approxi- 
mately a scalar multiple of Vj , where the scalar is 



'^u|x^/UjVj. yf course, we al'read) knew that behaved 
like a scalar multiple of the right eigenvector for 
large k, but now we have-a formula for ^he scalar. < 



fcxercise 6 . Return to the example of Section 2.3. Find an eigenvec- 
T 

tor of, A corresponding to a, = 1. Use this, together with the right 
eigenvector , 



v, - 2 



to find \ lm A*// * 



Let 



fl 6> 

a = ! r. * 

•a. Compute the dominant eigenvalue / of A, 

b. Find a positive eigenvector of A corresponding to 

T 

c. Find a positive eigenvector of A corresponding to 

d. find lin A k />^. 

k—i 1 

e. Describe the population distribution vector x, as k-«°. 
Exercise 8 . Carry out tj/e analysis of Exercise 7a-e for the matrix 

fo V s i] 



A = \h 0 0, 

10 5 * oj 

A^at will eventually nappen to a population whose growth is governed 
by such .15 matrix A* 




4. A HUMAN POPULATION EXAMPLE 



Although we-began this unit by referring to the U.S. 

population, all of the examples *so far have dealt v%ith 

fictitious peculations uhose matrices v%ere concocted to 

make the computations simple. Let's end the unit by 

returning to the U.S. population. Spec 1 f ical 1> , v*e will 

consider the population of U S. females in iy6". It is 

ecu en lent to diwde. the population into 10 age classes 

•J -4, 5-i 4 , Ij-14, 4 5 -4 j. Since the number of births 

to females over 50* is negligible v»e will consider only 

females aged 0-4J. The Leslie matrix consists largel> ,of 

zeros. The onl> nonzero entries are in the first rov. (the 

a 's, and m the sub-diagonal fthe b.'s) 1 These entries 
l & i 

are gr.en i; t laMc 2. 



TABIC 2 



Start Of 
\gc Interval 


I irst Rov* 


Sub 


-diagonal 

f V 


4 


ooooo 




.99694 


5 


.0010S 




99842 


lu 


08203 




99^85 


15' 


28849 


A. 


99071 


2(> 


. 3"S0 * 




99614< 


25 


.26478 




.99196 


3u 


. 14055 




9924 " 


55 


05857 




98875 










iu 


,01314 




98503 


13 


O0081 







Source Kevfitz and Flieger fsee References) 



I he don man t eigenvalue can be approximated using numeric a 
mathod^. It i^ approximately 1.0376. (Recall that each 
tine period la^ts 5 >ears f so /j = 1.0 3 "6 means that the 
total population *si:c l multiplied by approximately 



1 .0376^ every 5 years .) 
vector is ^ 



A corresponding positive eigen- 



\ 



'1 



/ 



f8.539 N * 
[8.204 
j :.394 
j - . 591 [ 
-.292i 

; ? .ooo I ' 

I6t-12! * • 
j 6 . 4 20 | 
'6.117! 
^5..95J 

Of course any positive scalar multiple of would do. The 
particular chofee^ of given here is useful because the 
entries represent percentages in the various age classes. 
For example, if 196" birth and survival rates persisted, 
then eventually approximately 8.5391 of the population 

'would be aged 0-4. 8.2041 would be aged 5-9, etc. 
Naturallv the percentages shown do not add up^ to \00*. 

'The remaining 281 [approximately) of the population would 
have age SO years or older,. It is interesting to note 
that*m 196" only about 241 of the population was SO year^ 
of age or older. It is well known that th>e average age 
of .the U.S. population is continuing to increase. 
Demographers believe that this increase will continue well 
into the twenty-first century, hence the concern about 
pension plans, etc. * 



Lxe r^ise 9 . 

Explain how to use high powers of a Leslie matrix A to approxi- 
mate left and right eigenvectors of A, (For a hint abou£ finding 
right* eigenvectors see the final sentence of Section 2 $.} * 

Based on your answer to part (a) show how to fvnd positive left 
and right eigenvectors for the Leslie matrix uf this section fU S. 
females, 1967 j using the computer program from Appendix A (or 
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* 4 

another program of your own writing). How could you use high powers of 
the Leslie matrix to find the dominant eigenvalue? 



5. RE I ERENCF^ , \j*> 



In thib unit we have given only a brief introduction 
to the Leslie model f or *populat lOn growth. There ^are many 
interesting extensions bf the model and uses of the model 
in related applications. If you want to read more about*, 
this model and see some applications to harvesting, you 
mi gut begin by examining 

Rorfes, C. and-H. Anton. Applications of Linear 

(P * T~~ 

* • Mgebra . John Wiley and Sons, 19//. 

i ^>ee 'par-t ) ei|l ar ly Chapters 6 , 9 and 10.) 

lor a di.^ussion of population models, including the Leslie ♦ 
model and others, you find these helpful 

P i e 1 o u , I . C . \n IntrQduction to Mathematical Ecology . 
John hi ley and^,§ons, 1969.^ ^ 

Poflard, J. H. Mathematical Models' for the Growth of 
Human Populat ions . Cambridge University Press, 

4 book which considers many interesting questions about 
-human populations is 

Keyfitz, N . Appl led Mathemat ical Demography . John 
hiley and Sons 
\ good source of data on human populations is 

Keyfitz, N . and W. Flieger. Population: Facts and 
, Methods of Demography . W. H. Freeman and Co., 
1971. 

The study of nonnegativ'e matrices is fascinating and you 
may wish to pursue it further. if so, you may. look at 
Chapter 2 of ( y-^ 

-Varga, R . S. Matrix Iterative Analysis . Prentice- 

Hall, 1962. ' 20 



O 1 



in thiN unit ut- have generall) dealt with 'smal 1 -sized 
matrices. lor such* raatr ice* we could easily find the eigen- 
valuer and eigenvector* and then use these to compute 
powers of the matrix, 'lor large matrices hand computation 
i- not practical. \ number of numerical methods are a vail - 
anle for computing eigenvalue* and eigenvectors of large 
matrices. To *ec more about the^e methods consult 

} o\% I . An Introduction to Numerical Linear Algebra . 
Oxford Uni\eisit> Press, 190 5 . 
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ANSWERS TO EXJERCISES 



(ISM] 

630 ! , 

! 

1 500 1 



r 1650 >l ■ , fl400| 



I 75 °j' 
I 325j 



[1575. 
j 700 
.[ J12.S) 



There are dore individuals in Class I than in 'Class II (roughly 
twice as many) and more in Class II than in Class III (again, about 
twice as many) . If we begin with 

•100' 



= j zoo j 

l 40U j 


, we nave 


tor x i , 


x : , .... x 6 






'1000- 


^250] 

1 ! 


f550] 


^625 1 . 


400 ] 


r 587.5 1 
1 


i 5 °' » 


,500] , 


1 i 
' 1* 


i j 
275 |, 


512 5 , 


200 | 


i lOOj 


i 25) 


(isoj 


1 62,5j 


137. 5j 


156. 25 j 



If »ou try several other starting vectors, you can probabLy con- 
vince yourself that the trend mentioned abov£ seems to hold 
regardless of which starting v«*;tor is chosen. 

Individuals in Class I produce 1 offspring on the average, %ndi- 
v iduaj ^s in Cla ss II produce 3 offspring on the average. Of the 
individuals in Class I, one-fourth can be expected- to suiv'ive 
until the next tame period 

, ¥* - 3 

detf'I-Aj «/*"-?- 3/4 = f 4 -|)f/ + ^) t so the eigenvalues are ^ 
and -i . 



and v-. 



t-ij 



K k ^ 

\ I > j approaches 

II I I 1 
II -1J 10 Oj ii -li 



are corresponding eigenvectors 



f 3/4 3/2] 



[1/8 1/lj 



The per cents are determined by v. . 0/7 or about 85 
Class I, 1/7 or about 14 will be in Class II 



will be in 



If | j is the initial age distribution vector, x^ is approxi- 
Ibj 



mately 



(1.5)' 



'3/4 3/2 
,1/8 1/4 



(1.5) k (ia ♦ ±b) 



Lach time k increases by 1, every entry in the age distribution 
^ vector is multiplied approximately by 1.5, so the population j 



3a, 



grows by 50%. 



>- ai 


" a 2 


-a 


" b l 


X 


0 


0 


" b 2 


\ 



- X* - a,A - b,a~A - b.b-a„ = 0. 
1 12 1 2 i 



There is one variation in sign (remember that all the a *s and 

1 s are nonnegative J . It follows from Descartes' rule of signs 
that there is exactly one positive root. You can show by induc- 
tion that the characteristic polynomial fpr the n x n matrix is 



,n \ n- J u A n-2 . u ,n-3 
} - a^A ^-b ] a ? A -b^b.a^A 



-b.b^-'-b .a . 
1 2 n-1 n 



Descartes' pule of signs again implies that there is one positive 
roo t . y 

b, y = bjX/A^ (note that Aj is nonzero); z = b^./A ^ ^b^x/Aj . 
If we set x = 1, we see that 



4. 



bl b 2 /xJ 



is a positive eigenvector corresponding to A^. Substituting 
these compo"nen<£ into the equation a^x + a 2 y +• a^z = A^x, we have 
2 

a l + a 2^1 / '^l + a 3^1^2^1 = ^1" This J- s equivalent to 

<$> -* — * 



A 1 - 



b 1 a 2 A 1 



^1^2 a 3 = ® y W ^ 1C ^ * s J ust result of 



substituting A^ into the characteristic equation. 







40' 




'320 




160' 


'1280' 




640' 


1 5j 


* 


40j 




20 > 




16oj ' 


80 , 


* 


640, 



The distribution vectors alternate between ones witjjt ratio 16:1 
and those with ratio 1 1. Every two periods the total population 
size grows by a factor of 4. The eigenvalues are 2 and -2, wit*h 
corresponding eigenvectors 
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'4 




'-4' 




and 




.1. 







thus * 



A k /2 k 



4 -4 
1 1 



1 0 



0 (-l)M 



4 -4 
1 1 



-1 



* k k 

Here (A^A^) = (-1) , which does not tend to zertfy but rather 
produces the oscillating effect. 

det(Xl - A T ) = det(AI T - A T ) = det((AI - A) T ) i= det(AI - A). An 
example; for the matrix ' 

0 8] 
A = I 5 * °J 

an eigenvector corresponding to eigenvalue 2 is 

ij ' «r 

T 

You can easily check that th^s vector is not an eigenvector of A . 



6. 



is an eigenvector of A corresponding to A^. 

4 8 8' 



•T T 1 
The limiting matrix is v^u^/u^v^ = 



2 4 4 
1 2 2 



* 7a, Eigenvalues are 2 and -1. 

b. An eigenvector of A corresponding to 2 is 

T 

c. An* eigenvector of A corresponding to 2 is 

k„k . T, T 



d. The limiting matrix for A /A^ is v ^ u j/ u i v i = g 



6 18 
1 3 



e. x^ behaves like a multiple o£ 



Each time k increases by 1, 



the entries in the age distribution vector are multiplied by , 
approximately 2. 

8a. Eigenvalues are *$(-! 1 i).. 
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On 



An eigenvector of A corresponding *to h is 



An eigenvector of A corresponding to ? is 



2 
2 
1, 

'l 
1 

UJ 



t k,,k 



The limiting matrix for A /X* is 
. /- i 



.4 .4 .4 
■4 .A .4 

.2 .2 .2 



x, behaves like a multiple of 
k 



2 ! 

2 | ' 



The entries in the age dis- 



tribution vector are multiplied by approximatefy 4 each time k 
increases by 1. The population will eventually die out (prac- 
tically, if not mathematically). 

k T 
For large k, A is approximately a scalar multiple of v jU^. But 

in the matrix v uT the columns are all multiples 6f v. and the 

11 ' r 64 . 

rows are multiples of u.. - Hence, any column of A is, an approxi- 

64 f 

mation to a right eigenvector of A and any row of A is an 

approximation to a left eigenvector of A. Each time k increases* 

k . > 

by 1, the entries of A are multiplied by approximately a.. 

k+1 r 

Therefore, to find we can take any entry of A (e.g., the 

program of Appendix -A could easily be modified to compute A ) and 

k 64 

divide by the corresponding entry of A (A m our program). 



APPENDIX A 



- Here is a BASIC prografi which computes powers of a 
matrix A. A sample run is also shown. 



LIST 

POWERS * . 

5 DIM A ( 1 0 , 1 0 ) 

1 0 INPUT N * 

1*S^!AT INPUT A ( N , N ), . ' 

20 PRINT ^ 

25 PR I NT ' 1 MATRI X A 1 ' * 

30 PRINT 

3 5 MAX PRINT A; 

40 FOR I =^1 TO 6 ^ * 

4£ MAT B = A* A 

50 PRINT M MATRIX A 1 1 j 2* * I 

55 PRINT ( ^ 

6 0 MAT 'PRINT ET| 

65 MAT A = B t ■ • „ 1 

70 NEXT I . 
75 END 

Ready * « 

run * *" 

POWERS 

? 0,1, 2, .5, 0.0,0-, ,5,0 ' . T 

MATRIX A 

-• * 

0 1 2 

* ' m » m 

a . 5 0 0 * # 

0.50 

Matrix a 2 
.5 1. 0 - * 

0,51 • 
• 2 5 0 a 

MATR f X A 4 * 

f 

.25 1 1 ' 

) 

^25 .,25 .5 * 
.125 .25 0 
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ERIC 



J 



MATRIX A 8 
.4375 .75 .75 
. 1 875 . 4375 . 375 
. 09375 . .1875 .25 

MATRIX A 16 
.402344 .796875 .796875 
.199219 . 40234 4 ^ .398438 
.996094E-1 .199219 .203125 

MATR^ A 32 
.400009 .799988 .799988 
.199997 .400009 .399994 
.$9*9985E-1 .199997 .200012 

MATRIX A 64 

.4 v 8 .8 

.2.4.4 
i 

' . 1 .2 .2 
ftEADY 
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APPENDIX B 



he will look once more at powers of a Leslie matrix A 

and show hov, the limiting behavior can be described using 

left and right eigenvectors. We WiTl continue to make 

those assumptions about A which are stated in SeLtion 3.4. 

7 

As before, let u, be a positive eigenvector of A corres- 

T 

ponding to the dominant eigenvalue A^, i.e., A u^ = '^^J 
alsc, let Vj be a positive eigenvector of A corresponding 

Let P be the square matrix v>hose 



»i . e . , Av , 



columns are made up of the eigenvectors of A, with 
being the first column,' i.e., 



rv r v ; 

0 , . . 

1 ... 



3 



• v n J 



-1 -if : 
p A p = p .v^v. 



P V 



o" 



i 0 j 



Also since P^AP- , * . j 

and since a diagonal matrix is its ov>n transpose, 



P X AP = (P ^P) 1 - PV(P j 1 



We will assunc, as before, that A is diagonalizablc, the final 
result is true .without this assumption. 
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,ERIC 



r 



-11 _ _ _ 
This neans that fP j = fu,u ? ...u ) where the u are 
y i *. n 1 

eigenvectors of A . u^ need not be the same as , but it 
will j?e a multiple of u ^ , say u^ = au^. 



It follows that P^UjU^.-.u j = P^P'^) 1 



T * ' T— 
Therefore P u^ = j.P u^ 



SO r UA = 



0 , 



fl 0 

I 

1 0 1 
I o * 6 

!?. 

I • t 
,0j 



a n^L t ir.pl les that u^P = f'i0...0j. 

Putting fB-1, and (B-2) together, we have 



0^ 



0 I 



Ij 



(B-2) 



1 T 

p W = 



T n -1 



'P f-ji. . .0, 

V 

' 0 i 



and u^Vj = u^PP = (a0...0j. ' 

10; 



fl IT*...' 



'1 01 
-0 0) 



and therefore, 



u i l 



= P' ip 

10 oj 



- 1 



Thus vve have te\ 



r 
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:;3 




and 



'nj/- 



0 



0 i 



0 _ 

.k' 



0 * 0 



-1 



0 1 

1 1 



Vl 
~ T 



2 / 1 



k / k i 



b 1 



o . p -i t 



n' lj 



As k-«, v*e see that A*/>j approaches v^Uj/u^v^ whi-ch is *hat 
we wanted to see. 
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